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ABSTRACT 

Bars in gas-rich spiral galaxies are short-lived. They drive gas inflows through their 
gravity torques, and at the same time self-regulate their strength. Their robustness 
has been subject of debate, since it was thought that only the resulting central mass 
concentrations (CMCs) were weakening bars, and only relatively rare massive CMCs 
were able to completely destroy them. Through numerical simulations including gas 
dynamics, we find that with the gas parameters of normal spiral galaxies, the CMC is 
not sufficient to fully dissolve the bar. But another overlooked mechanism, the transfer 
of angular momentum from the infalling gas to the stellar bar, can also strongly 
weaken the bar. In addition, we show that gravity torques are correctly reproduced in 
simulations, and conclude that bars are transient features, with life-time of 1-2 Gyr in 
typical Sb-Sc galaxies, because of the combined effects of CMCs and gravity torques, 
while most existing works had focussed on the CMC effects alone. 

Key words: galaxies: evolution - galaxies: spiral - galaxies: structure - galaxies: 
kinematics and dynamics. 
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1 INTRODUCTION 

In most barred galaxies, gas is concentrated on the lead- 
ing s ide of the bar (e.g., Ide Vaucouleurs fe de Vaucouleursl 
1963) , which is rep roduced in numerical simulations (e.g., 
lAthanassoulal FL992) . Then, the bar gravity torques make 
gas lose angular momentum, which initiates a gas inflow 
and fuels a central mass concentration (CMC). The bar is 
weakened by the CMC growth, because of escaping orbits 
(Hasa n fc Normanll9 90: Pfcn niger fc Norman 1990). Friedli 
(1993), iBerentzen et alJ il998t) and iHozumi fc Hernauisd 
(1999) report the dissolution of galactic bars with CMCs 
of mass 0.5 to 2 % of the disk mass. To explain the large 
fraction of barred g alaxies observed b oth at z ~ (e.g ., 
lEskridee eta!] 120021) and at z > 0.7 dSheth et al J 120031) . 
we have proposed that ba rs are dissolved and reformed 
dBournaud fc CombeJ2002l hereafter paper I) : our N-body 
simulations showed that bars are dissolved in 1-4 Gyrs in 
most galaxies, which we had attributed to the CMC growth. 
Yet, whether bars are really transient is still debated. Reso- 
nant rings are often observed in barred galaxies, which im- 
plies that gravitat ional torques are muc h larger than vis- 
cous torques (e.g.. lButa fc Combeslll99r]) . For a long time, 
gaseous simulations were too limited in resolution, induced 
too large viscous torques, so that such rings did not form. 
If viscous torques are still over-estimated, this may induce 
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unrealistic inflows of gas, so that CMCs are more massive 
and/or fueled more rapi dly, and the life-time of bars may be 
much under-estimated (IReean fc Teubenl 120041) . Moreover, 
even if the CMC fueling is realistic, the effects of the CMC 
may not be enough to fully dissolve the bar: recent sim- 
ulatio ns with more spatial resolution by IShen fc Sellwoodl 
(120041) - that do not include the whole gas response - have 
shown that bars are more robust against the growth of 
CMCs than what was believed before. 



Thus, the limited lifetime of bars found in paper I and 
other works could be an artifact of viscous torques or an 
over-estimation of the effects of the CMC. In this Letter, 
we show that the gas inflow in our simulations is actually 
initiated by gravitational torques, that are much larger than 
viscous torques. While most existing works had focused on 
the effects of the CMC growth, we show that another phe- 
nomenon can also lead to the destruction of galactic bars 
in a few Gyrs: the torques between the stellar bar and the 
gaseous arms also disturb the bar and can dissolve it. We 
then find that galactic bars are transient, even when the gas 
infall and the effects of the CMC are correctly treated and 
do not destroy the bar themselves: the growth of the CMC 
is not the only factor of bar dissolution. 



2 F. Bournaud, F. Combes and B. Semelin 



2 NUMERICAL SIMULATIONS 

We employ the N-body FFT code described in paper I, in- 
cluding star formation. The gravitational potential is com- 
puted on a Cartesian grid of size 512x512x64. The softening 
length and cell size are 75 pc, and the number of particles 
is 10 6 for each component (gas and stars). 

The dissipative dynamics of the ISM is modeled by the 
3-D sticky-particles code described in paper I (we use an 
elasticity parameter for cloud-cloud collisions @ r = fit — 0.75 
in this paper). Star formation and stellar mass-loss are also 
included. We use a visible (stellar and gaseous) mass of 
2 x 10 11 Mq. The disk truncation radius is 14 kpc and its 
scale-height is 1.2 kpc. The bulge and dark halo are Plum- 
mer spheres of radial scale-lengths 2 and 36 kpc respectively. 
The bulge-to-disk mass ratio is 0.25, and the dark-to- visible 
ratio inside the disk radius is 0.6. The initial gas mass is 
14.5 x 10 9 Mq (7.25 % of the visible mass). At t = 1.5 
Gyr, during the bar dissolution, the remaining mass of gas is 
11.3 x 10 9 Mq (5.5 % of the visible mass). We also run a con- 
trol simulation without gas. We compute the bar strength 
Pi as the maximum over radius of the m = 2 angular com- 
ponent of the gravitational torques, following the definition 
given in paper I, where this parameter was denoted 52 (see 
eq. 27 in paper I). 



3 MECHANISMS FOR BAR DISSOLUTION 

3.1 Gas inflows and bar dissolution 

The response of gas to the stellar bar is shown in Fig. 
Inside the corotation radius (CR, 6kpc) , gas is concentrated 
on the leading side of the bar. Then, as shown in Fig. [5] 
the bar gravity torques make gas lose 10-15% of its angu- 
lar momentum over one rotation, while the shear viscosity 
torques are much smaller: the gas inflow in our simulations 
is really initiated by the bar gravity torques, contrary to the 
hypothesis raised bv lReean fc Teubenl i2004h . 

The bar-driven gas inflow fuels a central mass concen- 
tration (CMC): we show in Fig. [^1 the increase in the mass 
inside radii 100 and 250 pc. At t = 2 Gyr, the mass radial 
profile shows a central peak of stars and gas that can be 
fitted by a Plummer CMC of mas s 2.8 x 10 9 M^ and radia l 
scale- length 90 pc. According to IShen fc Sellwoodl i2004l) . 
such a CMC cannot fully dissolve the bar. Yet, while the 
gravitational softening lengt h and the number of par ticles 
we used are close to that of IShen fc Sellwoodl i2004T) . the 
bar is here fully dissolved in less than 2 Gyrs: the residual 
values of P2 — 0.15 correspond to a weak oval distortion and 
spiral arms, but no bar is present. 

A surprising point shown in Fig. [3] is that the CMC 
growth occurs after the bar is already significantly weak- 
ened. The increase of the mass included inside radius 250 pc 
occurs between t = 1.2 and t = 1.8 Gyr. The bar strength 
is then smaller than 40% of its maximal value. As for the 
increase in the mass inside radius 100 pc, it occurs at the 
very end of the bar dissolution, from t = 1.4 tot — 1.8 Gyrs. 
At the opposite, the bar weakening is already significant at 
t = 1 Gyr: at this time, no strong increase in the mass has 
occurred, excepted at large scales of more than lkpc, but 
mass concentrations at such scales are totally inefficient to 
dissolve bars. It thus seems that the CMC growth is not the 
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Figure 1. Face-on plot of the gas density at t = 800 Myr. The 
dashed line represents the bar major axis, defined as the major 
axis of the most eccentric ellipse in an ellipse-fitting model of the 
stellar density. Gas is concentrated on the leading side of the bar 
inside the CR radius (6 kpc), with a shift of several degrees. 




0.0 0.5 1.0 1.5 2.0 

radius 



Figure 2. Angular momentum lost by gas in one rotation, in units 
of the inital angular momentum L(r), as a function of radius (in 
units of the CR radius at t = 900 Myr). The curves labeled "G" 
show the effects of the gravitational torques; they are mean values 
over periods of 100 Myr around t = 800 (dotted line), 900 (solid 
line), and 1000 Myr (dashed line). The curve labeled "V" shows 
the effects of viscous torques, that are negligible with respect 
to gravity torques. Inside the CR radius, gas loses 10—15 % of 
its momentum in one period. Torques are positive between the 
CR and the outer Lindblad resonance (OLR). Because of the bar 
slowing-down, the CR and OLR radii increase with time. The thin 
solid line shows the comparison with an SPH simulation with the 
same initial conditions at t = 900 Myr. Gravity torques are of the 
same order with this different code. 

only factor that is responsible for the destruction of the bar. 
Furthermore, the bar dissolution is not an artifact of a too 
long time step: we have divided the time step by factor 8, 
and found a similar result (see Fig. . 

3.2 Effects of the CMC 

To confirm that the CMC growth is not the only phe- 
nomenon responsible for the bar dissolution, we have run 
a purely stellar simulation with the same initial conditions, 
and ad ded an analytical CMC. T his model is similar to that 
used bv lShen fc Sellwoodl <2004 . The final CMC is a Plum- 
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Figure 3. Evolution of the bar strength in the main simulation 
(solid line labeled "P2") and in the simulation with an eigth times 
shorter time step (dashed), and increase in the central masses 
inside radii 250pc (dashed line) and 100 pc (solid line). The ini- 
tial masses have been subtracted, so that these curves only show 
the bar-induced growth of a CMC. All the quantities have been 
rescaled to 1 at their maximum value. 
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Figure 4. Evolution of the bar strength in several models: "gas" : 
complete simulation including gas - "control": stellar simulation 
without gas — " CMC" : stellar simulation without gas and addition 
of a growing CMC - "gas-CMC": complete simulation including 
gas and artificial removal of the CMC growth - "torques": stellar 
simulation without gas and addition of gravity torques represent- 
ing those exerted by gas on the stellar bar. 



mer sphere corresponding to the CMC we had fitted in the 
complete simulation (see previous Section). Its radial scale- 
length is 90 pc, and its mass is grown up to 2.8 x 10 9 Mq 
following the variations of M(100pc) shown in Fig.|3] In this 
model, the bar is not fully dissolved, but only partly weak- 
en ed (see F ig. |H. T his is i n agreement with the conclusion 
of lShen fc Sellwoodl lj2004T) that the bar can survive to the 
growth of such a CMC. 

We have also run a simulation in which we include 
the gas response, but suppress artificially the growth of the 
CMC: every gas particle that is found below radius 100 pc, 
with a velocity inferior to the circular velocity at this radius, 
is artificially suppressed, unless already present at t = 0. We 
do not suppress particles that are just crossing the central 
100 pc with a large velocity, but only those that have defini- 
tively fallen in this region. As shown in Fig. [I] the bar is 
dissolved, and the evolution is very similar to that of the 
complete simulation. This result definitively proves that a 
phenomenon different from the CMC growth intervenes in 
the dissolution of the bar. This phenomenon is related to 
the gas response, since it is not observed in the purely stel- 
lar simulation (Fig.^J. 

3.3 Effects of gravity torques 

Inside the CR, the bar has been shown to exert important 
gravitational torques on the ISM. In the same time, the lead- 
ing arms of gas exert torques on the stellar bar. The mean 
torque undergone by stars inside the bar is positive, but the 
real torques vary with radius and azimuth. This could be the 
reason why the bar is dissolved in the simulation where the 
CMC has been artificially suppressed: gravitational torques 
between the stellar bar and gaseous arms are still present in 
this model. 

To confirm the role of gravity torques, we have run a 
purely stellar simulation in which gravitational torques ex- 
erted by gas have been artificially added. To this aim, in 
the complete simulation, we compute the amplitude Ai and 
phases 4>i of the tangential forces _F,| as exerted by gas in the 
Fourier decomposition 



ff(r,0) =^^(r)sin(i(0-^)) (1) 

i 

where 6 is computed in the frame of the bar. We record the 
mean values < Ai > and < cj>i > between 0.8 and 1 Gyr. 
Then, in a purely stellar simulation, we analytically add the 
gravity torques given by 

F| as M)= < A i > W sin (*(*- < h >)) ( 2 ) 

1=2,4,6,8 

where 9 in still computed in the frame of the stellar bar. 
The amplitude of these torques is realistic, since < Ai > 
and < 4>i > are given by the complete simulation, even if 
they do not follow any temporal evolution. We apply no 
torques before 800 Myr and increase them progressively up 
to their final amplitude between 800 and 900 Myr. In this 
model, that only differs from the purely stellar simulation by 
the addition of these gravitational torques, the bar is fully 
dissolved (see Fig. , an d its dissolution is nearly as rapid 
as in the complete simulation. 

Then, the gravity torques exerted by infalling gas on 
the bar can strongly weaken it. In this first simulation, this 
process dissolves it, while the role of the CMC is negligible. 
We have run another simulation with an initial gaseous mass 
equal to 4.5 % of the visible mass, and a bulge-to-disk mass 
ratio of 0.33, and made the same models with the CMC alone 
and the torques alone. The results are shown in Fig. |S] In 
the complete simulation, the bar is fully dissolved, while the 
growth of the CMC only partly weakens it. Here again, the 
gravitational torques play a major role in the bar dissolution. 
They are not strong enough to fully dissolve the bar alone, 
but the combined effects of the CMC lead to a complete 
dissolution. This second case confirms the role of gravity 
torques in the dissolution of bars, even is the CMC plays a 
role, too. 



4 DISCUSSION AND CONCLUSION 

In our simulations, viscous torques are negligible with re- 
spect to gravitational torques and do not directly initiate 
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Figure 5. Same as Fig. HI with an initial gas-to-visible mass 
ratio of 4.5% and a bulge-to-disk mass ratio of 0.33. The bar is 
dissolved in the complete simulation (curve labeled "gas"). The 
effect of gravity torques (curve labeled "gas-CMC") participate 
much to the bar weakening, but the CMC also plays a role. 



the gas inflow. Yet, the phase shift between the bar and 
gas, and the magnitude of the result ing gravity torques, are 
large ly dependent on the viscosity llvan Albada fc Roberts! 
1981). Then, whether gravity torques in our simulations are 
realistic or not has to be checked. 

First, we have checked that our results are not depen- 
dent on the box size used to compute particles collisions 
in the sticky-particles scheme: this parameter controls the 
distance along which particles interact, influences the nu- 
merical viscosity, thus the gravity torques may depend on 
it. For box sizes ranging from 10 pc to 350 pc, we find that 
changes in this parameter does not lead to large changes in 
the gravity torques, and the bar life-time is not dependent 
on this parameter, as shown on Fig.|S] It may seem surpris- 
ing that changing this box size parameter does not influence 
the bar life-time at all. In fact, as shown by Fig- El the grav- 
ity torques are actually dependant on the sticky particles 
parameters (the box size, and the elasticity parameter j3 for 
cloud-cloud collisions) : gravity torques increase with the box 
size, and decreases with the elasticity parameter. But there 
is a range of parameters where gravity torques are nearly 
constant, for elasticity parameter f3 between 0.65 and 0.85, 
and box sizes between lOpc and 350 pc. For larger box sizes 
and/or smaller /3, gravity torques become larger, but the 
medium is very cold (with a typical sound speed smaller 
than 5 km s _1 ) and clumpy, instead of forming spiral arms: 
these situations are not realistic, so they are ruled out. At 
the opposite, large f3 or small box sizes lead to situations 
with smaller gravity torques, but the dissipation rate in the 
ISM is not realistic, the sound speed in the ISM is found 
to be larger than 10 or even 15 km s" 1 and the Toomre 
parameter for gas is larger than 1.5: as a result, gas only 
forms very smooth spiral arms that disappear after a few 
dynamical times. Thus, the only acceptable parameters are 
0.65 < P < 0.85 and box sizes from 10 to 350 pc (for our 
number of gas particles). In this range, there is an equilib- 
rium between dissipation in the ISM and kinematical heat- 
ing by density waves: for instance a smaller (3 leads to more 
dissipation, so the first spiral arms are stronger and heat the 
medium more. This explains why physical properties are not 
strongly affected by the gas dynamical parameters in this 
particular range. Then, provided that we choose parameters 




Figure 6. Evolution of the bar strengh with different box sizes 
in the sticky-particle scheme, for the first simulation detailed be- 
fore. The thick solid line represent this main simulation (see also 
Fig. 1), with a box size of 50 pc. Other lines represent box sizes 
from 10 pc to 350 pc. The bar dissolution is not largely affected by 
this parameter. The inset shows the gravity torques for collision 
box sizes of 10, 50 and 350 pc (see Fig. 2 for details). The torques 
are slightly larger when the box size is larger, but the change in 
gravity torques is rather small: it does not exceed the variations 
observed in 100 Myr (see Fig. 2). 
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Figure 7. Variations of the gravity torques exerted on the gas, 
as a function of the box size (solid curve, bottom axis, with an 
elasticity parameter (3 = 0.75), and as a function of the elastic- 
ity parameter (dashed curve, top axis, with a box size of 50 pc). 
the quantity AL/L, plotted in absolute value, is defined as in 
Sect. 3.1. Gravity torques are not strongly affected by these pa- 
rameters in the range 10— 400pc for the box size and j3 =0.65-0.85 
for the elasticity parameter. This explains why the bar evolution is 
not strongly affected by a change in these parameters, as shown in 
Fig. 6. Outside this particular range, torques are more affected by 
these two parameters, but the ISM is either unrealistically cold or 
unrealistically stable : only parameters in the range where gravity 
torques are not strongly affected are acceptable. 



resulting in a realistic model for the ISM (not too cold or 
too hot, i.e. a medium in which the main instabilities are 
spiral arms), the physical results are not strongly affected 
by the sticky particles parameters, so they are robust. 

As a second verification, we have estimated the gravity 
torques exerted by stellar bars on gaseous arms in a sam- 
ple of observed barred galaxies (see a forthcoming paper for 
details). We show on Fig.|S]the results for three well-known 
barred galaxies, NGC 1365, NGC 7479, and M 100: for these 
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Figure 8. Gravity torques exerted by the bar on the gaseous 
arms in NGC 1365, NGC 7479 and M 100, as function of the 
radius (in units of the CR radius), for several tracers of the ISM 
(NGC 1365: Ha and dust emission at 7[im - NGC 7479: CO and 
Ha - M 100: Ha and dust absorption traced by the B-V index). 
The gravitational potential of bars is estimated from I-band im- 
ages. The CR radius has been defined as the radius where the 
averaged torques become positive. Gas inside the CR loses 10 
to 20% of its momentum over one period: this is fully compati- 
ble with our simulations, which proves that the gravity torques 
between the bar and gas are realistic in our simulations. 



galaxies as for other barred ones that we have studied, inside 
the CR, gas loses typically 10-20% of its angular momentum 
in one rotation. This result is independant of the tracer of 
interstellar gas that we choose (see Fig. |HJ . This is compati- 
ble with our simulations (Fig. 2). Thus, the gravity torques 
are correctly reproduced in our sticky-particles simulations, 
that provide a realistic estimation of the life-time of bars. 

Finally, we also ran a simulation with the same ini- 
tial conditions as our firs t run, with the tree-SPH code by 
ISemelin fe Combesl (120021) . using 2 x 10 5 particles, isother- 
mal initial conditions at 10 4 K, and adiabatic evolution. As 
shown in Fig. 2, the gravity torques between the stellar bar 
and gas have the same order of value as in sticky-particle 
simulations, and in observations. 



While the dissolution of bars had always been at- 
tributed to the growth of CMCs, we find here that an- 
other process contributes to their destruction. Even if grow- 
ing CMCs do not com pletely dissolve most galactic bars 
iShen fe Sellwoodll2004l) . the gravitational torques exerted 
by gas on the bar can strongly weaken it: these torques 
are positive, which make orbits become rounder; they also 
induce a phase shift that depends on the radius, which 
misaligns stellar orbits, so that the bar is strongly weak- 
ened. In the two simulations with parameters typical of 
Sb-Sc galaxies, the cycle of bar dissolution takes about 2 
Gyrs (for the mass and radius of the Milky- Way) . The 
supply of angular momentum from the ISM to the bar, 
which dissolves it, occurs much faster than the removal 
of momentum by the darrk halo, t hat could trigger th e 
bar over timescales of several Gyrs (|AtJuirias i s i ouki|j2002j). 
The constant bar fraction up to z ~ 1 llSheth et al.l 2003; 
Ijogee et al.lll99ollElmegreen et alJl2005tl cannot then be in- 
ter preted in terms of r obust , long-lived bars. As suggested 
bv lSellwood fe Moord dl999h . the bar could be renewed ei- 
ther by galaxy interact ions, or by accretion of external gas. 
iBerentzen et alJ i2004) have shown that galaxy interactions 
can reform bar only in gas poor galaxies, i.e. not in most 
spirals. This indicates that gas accretion plays a major role 
in the high fraction of barred galaxies, which can accoun t in 
detail for observed properties of bars iBlock et alj|2002l) . 
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